Thermodynamic driving force of formation of coherent three-dimensional islands in 

Stranski-Krastanov growth 
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The formation of coherent three-dimensional islands in highly mismatched epitaxy is discussed 
in terms of the traditional concept of wetting. It is shown that the wetting layer and the 3D 
islands represent different phases which cannot be in equilibrium with each other. The transfer 
of matter from the stable wetting layer to the 3D islands is thermodynamically unfavored. The 
experimentally observed critical misfit for coherent 3D islanding to occur and the coexistence of 
pyramids with discrete heights of two, three, four... monolayers can be explained assuming that 
the thermodynamic driving force of formation of coherent 3D islands on the surface of the wetting 
layer of the same material is the reduced average adhesion of the islands to that layer and that the 
islands height is a discrete variable. 

PACS numbers: 68.55.Jk, 68.35.Md, 68.35.Np, 68.66.Hb 
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The growth of thin epitaxial films usually takes place 
far from equilibrium. Nevertheless, thermodynamic con- 
siderations are a necessary step for understanding of the 
process. Of particular interest is the thermodynamic 
driving force (TDF) which is responsible for one or an- 
other mechanism of growth. While this question is well 
understood in terms of wetting of the substrate by the 
overgrowth in the cases of island or Volmer- Weber (VW) 
growth .and, layer-by-layer or Frank- van der Merwe (FM) 
growth£ratl the Stranski-Krastanov (SK) growth (3D is- 
lands on top of a thin wetting layer) is far from being 
clarified. The reason is that the SK growth is in fact a 
growth of material A on the same material A, which ther- 
modynamically requires the formation and growth of 2D 
rather than 3D islands. This is particularly true, in the 
case of the coherent Stranski-Krastanov growth,cl where 
dislocation-free 3D islands-|are strained to the same de- 
gree as the wetting layer .aaa This is the reason why it is 
widely accepted that the energy of the interfacial bound- 
ary betwieen the 3D islands and the wetting layer is equal 
to zero.Q Although this energy is expected to be small 
compared with that of the free crystal faces ,□ it should 
not be neglected since this is equivalent to the assump- 
tion that the islands wet completely the wetting layer. 
The latter rules out 3D islanding from a thermodynamic 
point of view.EI 

The need for a thermodynamic analysis arises also from 
the experimental observations of a crit ical misfit for co- 
herent 3D islanding to occur and the simulta- 
neous presence of i slan ds of different thickness which vary 
by one monolayer .liSLa The existence of a critical misfit, 
as well as of stable two, three or four jaapnolayers thick 
islands, do not follow from the tradeoffllJ 
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between the cost of the additional surface energy and 



the gain of energy due to the elastic relaxation of the 3D 
islands relative to the wetting layer (V, 7 and £0 are the 
islands volume, the specific surface energy and the lattice 
misfit, respectively, and C and C" are constants). 

It was recently suggested that the TDF for coherent 
3D islanding is_the incomplete wetting of the substrate 
by the islands,E£l rather than the elastic relaxation of the 
material in the islands. The incomplete wetting is due to 
the displacements of the atoms near the island edges from 
the bottoms of the corresponding potential troughs pro- 
vided by the wetting layer. This results in a series of criti- 
cal volumes at which the monolayer high islands become 
unstable against the bilayer islands, the bilayer islands 
against the trilayer islands, etc. The misfit dependence of 
the first critical size N12 for the mono-bilayer transforma- 
tion displays a critical behavior in the sense that coherent 
3D islands can be formed at a misfit higher than some 
critical value. Below this value the film should grow in a 
layer-like mode until misfit dislocations are introduced to 
relieve the strain. However, the approximation used by 
the authors, whkh_is based on the ID model of Frenkel 
and KontorovaJ13ila was unable to describe correctly the 
individual behavior of atoms inside each layer, since it as- 
sumes a potential with a period given by the average of 
the separations of atoms (considered frozen) in the layer 
underneath. Although this model gives qualitatively rea- 
sonable results concerning the energy of the islands, it is 
inadequate to calculate, in particular, the average adhe- 
sion energy of the islands to the wetting layer. 

In the present report we recollect some simple thermo- 
dynamic aspects of the epitaxial morphology based on 
the traditional concept of wetting and consider the co- 
herent SK growth from this point of view. The same 
concepts were in fact advanced by Stranski in his model, 
admittedly very peculiar, of a monovalent ionic crystal 
K + A~ on the surface of an isomorphous bivalent crystal 
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Wc then support our thermodynamic consi- 
derations by numerical calculations making use of a sim- 
ple minimization procedure on the same atomistic model 
in 1 + 1 dimensions (length + height) as in Ref.©. 
The 3D islands are represented by linear chains of atoms 
stacked one upon the other,E3 the islands height being 
thus considered as a discrete variable which increases by 
unity from one. The latter is of crucial importance as the 
aspect ratio of the 3D islands is usually of the ordetxif 0.1 
and the height is of the order of 10 monolayers. EE The 
atoms interact through an anharmonic Morse potential 



V(x) = Vole-^'-Q - 2e- &{x ~ b) ] 
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The total interaction energy as well as its derivatives 
with respect to the atomic coordinates, i.e. the forces, 
are calculated. Relaxation is then performed iteratively 
by allowing the atoms to displace in the direction of the 
forces until these fall below some negligible cutoff value. 
We consider interactions in the first coordination sphere 
in order to mimic the directional bonds that are char- 
acteristic for semiconductors Bj Inclusion of further co- 
ordination spheres alters only minimally the numerical 
results. The substrate (the wetting layer) is assumed to 
be rigid. 




Number of monolayers 



FIG. 1: Schematic dependence of the film chemical pot- 
ential on the film thickness in number of monolayers for the 
three modes of growth: VW - Volmer- Weber, SK - Stranski- 
Krastanov, and FM - Frank-van der Merwe. The dashed line 
gives the chemical potential of the unstable wetting layer. 

A mother phase (a vapor) and a new phase (e.g. a 
strained planar film or unstrained 3D crystals) are in 
equilibrium with each other when their chemical poten- 
tials are equal. A transition from one phase (mother or 
new) to another takes place when the chemical poten- 
tial of one of the phases becomes smaller than that of 
the other. The TDF for this transition is the difference 
of the chemical potentials of both phases at the given 
pressure and temperature. The TDF which determines 
the occurrence of one or another mechanism of epitaxial 
growth (growth from vapor of a strained 2D layer or 3D 
islands) is the difference A/j, — /i(n) — /Zgj-, of the chem- 
ical potential, /x(n), of the overlayer which depends on 
the film thickness measured in number n of monolayers 
counted from the interface, and the chemical potential, 



/ig D , of the bulk 3D crystal of the same materially The 
thickness dependence of fi(n) originates from the thick- 
ness distribution of the misfit strain and, on the other 
hand, from the interaction between the deposit and the 
substrate, which rapidly decreases with the distance from 
the interface (-Eab —* i?AA)-a'EiJ 

If we deposit a crystal A on the surface of a crystal B, 
A/j, can be written in terms of the interatomic energies 
per atom, Eaa and Eab, required to disjoin a half-crystal 
A from a like half-crystal A and from an unlike half- 
crystal B, respectively :E3 



(i(ri) = m§ d + [E AA - E A B(n)} 
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The adhesion energy Eab includes in itself the thick- 
ness distribution of the strain energy and the attenuation 
of the bonding with the substrate. <2> = 1 — Eab/ Eaa is 
the adhesion parameter which accounts for the wetting 
of the substrate by the overgrowth. Ea_ &) is equivalent 
to the familiar 3-er criterion of Bauer .ErEM 

As follows from (^) the parameter 4> = A^l/Eaa is 
equal to the TDF for occurrence of one or another mode 
of growth relative to the cohesive energy Eaa ■ In the 
two limiting cases of VW (0 < <P < 1) and FM growth 
(<£ < 0, £o ~ 0), Afj, tends asymptotically with increas- 
ing film thickness to zero from above and from below, re- 
spectively, but changes its sign in the jcase of SK growth 
(<£ < 0, £o 0), as shown in Fig. |l]B'EJ Consider now 
Fig. ^ in terms of equilibrium vapor pressures instead of 
chemical potentials. Although the connection is straight- 
forward, (/i oc InP), such_a consideration gives a deeper 
insight into the problem. FJ Thus, as long as fi(n) < /ig D 
a thin planar film can be deposited at a vapor pressure P 
that is smaller than the equilibrium vapor pressure, Pq, 
of the bulk crystal, but is larger than the equilibrium va- 
por pressure Pi of the first monolayer, i.e. Pi < P < Pq. 
In other words, a planar film can be deposited at under- 
saturation Afi — kTln(P/Po) with respect to the bulk 
crystal. The formation of 3D islands {ii{n) > fJ,^) re- 
quires P > Pq, or a supersaturation with respect to the 
bulk crystal. 

Applying the above considerations to the SK growth 
leads unavoidably to the conclusion that the 3D is- 
lands and the wetting layer represent necessarily different 
phases and thus have different chemical potentials. The 
reason is that the two phases are in equilibrium with 
the mother phase (the vapor) under different conditions 
which never overlap. The wetting layer can be in equilib- 
rium only with an undersaturated vapor phase (P < P ), 
while small 3D islands can be in equilibrium only with a 
supersaturated vapor phase (P > Pq). The dividing line 
is thus A/i = fcTln(P/Po) = at which the wetting layer 
cannot grow thicker and the 3D islands cannot nucleate 
and grow. Hence the wetting layer and the 3D islands 
can never be in equilibrium with each other. 

It follows from the above that the derivative, dAE/dV, 
of the energy of the 3D islands relative to that of the wet- 
ting layer, gives the difference of the chemical potentials 
of the wetting layer (the dashed line in Fig. P and the 
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chemical potential of the 3D islands. In other words, it 
represents the difference of the supersaturations of the 
vapor phase with respect to the wetting layer and the 
3D islands. As the thickness and the energy of the wet- 
ting layer depend on the misfit it would be more suit- 
able if one chooses as ar*eference the bulk crystal rather 
than the wetting layer.E3 Transfer of material from the 
stable wetting layer (/iwL < A*3d) t° the 3D islands is 
connected with increase of the free energy of the system, 
and therefore, is thermodynamically unfavored. A pla- 
nar film thicker than the stable wetting layer is unstable 
and the excess of the material can be transferred to the 
3D islands if the necessary thermal activation exists. 




G9 eef 



15 20 25 30 

Atom Nr. 



FIG. 2: Vertical displacements of the atoms of the base chain 
of a coherent (a) and a dislocated (b), 3 monolayers-thick 
island, given in units of the lattice parameter of the wetting 
layer and measured from the bottoms of the potential troughs 
provided by the homogeneously strained wetting layer. The 
misfit £o amounts to 7% and the islands contain 30 and 34 
atoms in their base chains, respectively. 

We focus our attention on the adhesion between the 
3D islands and the wetting layer. Fig. ||is an illustration 
of the difference (and resemblance) between the classi- 
cal and the coherent SK mode. In the calculations, the 
sizes of the base chain of the island (30 and 34 atoms, 
respectively) have been chosen just below and above the 
critical size for introduction of misfit dislocations at the 
given misfit of 7%. As seen, in both cases the 3D is- 
lands loose contact with the wetting layer. The vertical 
displacements are largest at the chain's ends in the co- 
herent SK mode and around the dislocation cores in the 
classical case, but the physics is essentially the same. The 
mean adhesion parameter increases with the islands' 
height and saturates beyond several monolayers (Fig. |). 
In our model, # is calculated as the adhesion energy be- 
tween island and wetting layer at the given misfit (7%) 
minus the corresponding value for zero misfit. It can 
be seen that the compressed overlayers exhibit a greater 
tendency to coherent SK growth than expanded ones as 
expected, due to the anharmonicity of the potential (^|). 

The same tendency is seen more clearly in Fig. ^ which 
shows in fact the dependence of the TDF for formation of 
coherent 3D islands, A[i, on the lattice misfit. The lat- 
ter remains close to zero for expanded overlayers but in- 
creases steeply beyond approximately 5% in compressed 
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FIG. 3: Mean adhesion parameter <2> as defined by Eq. (1) 
as a function of the islands' height in number of monolayers 
for positive and negative values of the misfit of absolute value 
of 7%. Coherent islands of 14 atoms in the base chain were 
considered in the calculations. 



overlayers. This behavior agrees well with the misfit de- 
pendence of the critical size, N\2, for the mono-bilayer 
transformation to occur, as shown in Fig. |5[ where a 
steep rise of N12 with decreasing absolute value of the 
misfit is observed only in compressed overlayers. Note 
that it is-less sharp for expanded ones opposite to earlier 
finding.N 
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FIG. 4: Mean adhesion parameter of one monolayer high, 
coherent islands as a function of the lattice misfit. The islands 
contain 20 atoms. Data for both positive and negative misfits 
are shown in one quadrant for easier comparison. 

We discuss now the discrete character of the height 
of the 3D islands. The experimentally observed volume 
of the auantum dots varies roughly from 20000 to 50000 
atoms.Qty Typical values of the aspect ratio ottbc islands 
height and half-base are of the order of O.l.B'E-3 Thus a 
pyramid with a base edge of 100 atoms and aspect ratio 
0.1, and containing 22000 atoms, is only 5 monolayers 
high. The addition of 14400 atoms (a new base plane of 
120 x 120 atoms) requires only one more atomic plane. 
Calculations of the energy of islands having a shape of 
a frustum of a pyramid are usually performed assum- 
ing implicitly that the lengths of the lower, R, and the 
upper, R' , bases, and in particular the height h, are con- 
tinuous variables. Eq. (|l]) is obtained by using the Ter- 
soff approximation which neglects the gradient of strain 
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in a direction normal to the surface plane together with 
h <C Ru, and assuming h 3> c where c is the atomic 
distances This would be correct if the crystals contain 
at least several million of atoms. 
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FIG. 5: Misfit dependence of the critical size N12 (in number 
of atoms) for positive and negative values of the lattice misfit. 
The curves are shown in one quadrant for easier comparison. 

We conclude that the wetting layer and the 3D islands 
represent different phases which cannot be in equilib- 
rium with each other, and the SK morphology is a re- 
sult of the replacement of one first order phase transition 
(vapor - wetting layer) by another first order transition 
(vapor - 3D islands). The transfer of matter from the 



stable wetting layer to the 3D islands is thermodynam- 
ically unfavored. The experimental observations menti- 
oned above can be explained on the base of two assump- 
tions: the thermodynamic driving force for the coherent 
3D islanding is the incomplete wetting and the height 
of the 3D islands is a discrete variable varying by one 
monolayer. This leads to the results that (i) monolayer 
high islands with a critical size appear as necessary pre- 
cursors for 3D islands, (ii) the 2D-3D transition takes 
place through a series of intermediate states with dis- 
cretely increasing thickness that are stable in separate 
intervals of volume (see Refs. ( |T4] , p^| )) (iii) there ex- 
ists a critical misfit below which coherent 3D islands are 
thermodynamically unfavored (see Refs. ( |l^ , |lT| , ^2|Jl^ ) ) 
and the misfit is accommodated by misfit dislocations 
at a later stage of the growth. Compressed overlay- 
ers show a greater tendency to 3D clustering than ex-. 
panded ones, in agreement with experimental results.£2l 
Result (i) explains readily why the volume distribution of 
InAs/GaAs self- assembled quantum dots agrees well with 
the scaling functionsJor two-dimensional submonolayer 
homoepitaxy model. cB 
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